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Abstract 

We present a theory of spin, electronic and transport properties of a few-electron lateral trian- 
gular triple quantum dot molecule in a magnetic field. Our theory is based on a generalization of 
a Hubbard model and the Linear Combination of Harmonic Orbitals combined with Configuration 
Interaction method (LCHO-CI) for arbitrary magnetic fields. The few-particle spectra obtained as 
a function of the magnetic field exhibit Aharonov-Bohm oscillations. As a result, by changing the 
magnetic field it is possible to engineer the degeneracies of single-particle levels, and thus control 
the total spin of the many-electron system. For the triple dot with two and four electrons we find 
oscillations of total spin due to the singlet-triplet transitions occurring periodically in the mag- 
netic field. In the three-electron system we find a transition from a magnetically frustrated to the 
spin-polarized state. We discuss the impact of these phase transitions on the addition spectrum 
and the spin blockade of the lateral triple quantum dot molecule. 

PACS numbers: 73.21.La,73.23.Hk 
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I. INTRODUCTION 



There is currently interest in developing the ability to control and manipulate the total 
spin of individually localized interacting electrons as a prerequisite for solid-state nanospin- 
tronic and quantum information applications.-*^'^'^'^ Precise control over the number and 
spatial location of carriers can be achieved by confining them in lateral gated quantum 
dot devices, of which the single, ^'^ double,-'^'^'^^'^^ and triple^i^^ii^ii^ quantum dots have 
already been demonstrated. In particular, Gaudreau et alM^ reported a controlled charg- 
ing of a lateral triple quantum dot (TQD) molecule with = 1 — 6 electrons, with the 
ability to control the population of each dot independently. Preliminary experiments on 
quantum dot molecules in external magnetic field by Gaudreau et alM and by Ihn et ai— 
showed signatures of Aharonov-Bohm (AB) oscillations, indicating coherent coupling be- 
tween the constituent dots. In this work we present a theory of the magnetic field effect 
on the electronic, spin, and transport properties of an isolated triple quantum dot molecule 
with controlled number of electrons A^ = 1 — 6. 

Previous theoretical descriptions of isolated lateral multi-quantum dot devices in a mag- 
netic field focused on quantum dot molecules with one electron per dot using Hubbard, exact 
numerical diagonalization, and spin Heisenberg model.— i2S They showed magnetic field in- 
duced corrections to the Heisenberg model due to chiral spin interactions. Furthermore, 
for three dots in a triangular structure, one electron each, they established magnetic fielc- 
induced transitions from a lowest-energy spin doublet with total spin S = 1/2 (Ref. l2ll ) 
to S* = 3/2 spin polarized state.— 122 Spin transitions in isolated lateral multi-quantum dot 
devices with large electrons numbers have also been studied using spin density functional 
theory by Stopa et alM There has also been significant interest in triple quantum dots in 
triangular configuration connected to the leads. Using the Hubbard model the effects of the 
magnetic field on the conductance through an empty and singly occupied triple dot were 
studied, with focus on the interplay between the Kondo physics, symmetries, and the AB 
osciUations.^^^^^ 

The aim of this work is to study the magnetic field dependence of the electronic properties 
of the lowest electronic shell of a triangular triple quantum dot molecule filled with A^ = 
1 — 6 electrons, extending in this way our previous work^'' to finite magnetic fields. This is 
accomplished by both the analysis of the Hubbard model and by the development of a new 
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computational tool. The new microscopic tool combines (i) a calculation of single particle 
states as a linear combination of harmonic orbitals (LCHO) localized on each dot, with a 
proper gauge transformation allowing for convergent results as a function of the ratio of 
the magnetic length to the inter-dot separation, with (ii) configuration-interaction approach 
(CI) to the many-electron problem. These techniques have allowed us to analyze the spin 
and electronic properties as a function of the magnetic field and the number of confined 
electrons (A^ = 1 up to 6). We derive the magnetic-field evolution of the one-electron 
spectrum and show the existence of degeneracies at multiples of half fiux quanta threading 



the area of the TQD, in agreement with Ref. l23|. The magnetic field-engineered degeneracies 



of single-particle levels, combined with electron-electron exchange and correlations, allow for 
the control of the total spin of the many-electron system. For example, we show total spin 
oscillations due to the singlet-triplet transitions occurring periodically in the magnetic field 
for two and four electron molecules. In the three-electron system we find the magnetic field- 
induced transition from a magnetically frustrated to the spin-polarized state, in agreement 



with Refs. 
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|20|. We discuss the impact of these spin transitions on the addition spectrum 
as measured using charge spectroscopy, and predict the appearance of spin blockade in the 
transport through TQD molecule. 

The paper is organized as follows. In Sec. [Ill we present details of the Hubbard and 
LCHO-CI approaches. In Sec. IIIII we calculate the electronic structure of the triple dot 
filled with = 1 to 6 electrons as a function of the magnetic field. The discussion of the 
charging diagram and addition amplitudes is presented in Sec. HVl The paper is summarized 
in Sec. |Vl 

II. THE MODEL 

A schematic picture of the TQD studied in this work is shown in Fig. [TJ^a). This system 
is an approximation of the lateral gated TQD device, in which the three potential minima 
are created electrostatically by metallic gates. Such device has been studied theoretically 



in Refs. 
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and is related to the system demonstrated experimentally by Gaudreau et 
alr^ Figure [11(b) shows the TQD electrostatic potential generated by a model arrangement 
of gates enclosing the area of the device (not shown) together with additional gates (shown 
as white regions) used to establish the potential barriers between the dots. By selective 
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tuning of the voltages it is possible to bring the three dots into resonance, i.e., match the 
energies of the lowest single-particle orbital of each potential minimum. The resulting TQD 
molecule can be then filled controUably with N electrons, starting at = 1, in the presence 
of a magnetic field B = [0, 0, B] applied in the direction perpendicular to the plane of the 
system. 

A. Hubbard model 

We have shown previouslyi^ that the electronic properties of the molecule with few con- 
fined electrons (A^ = 1 to 6) can be understood in the frame of the Hubbard model. Assuming 
one orbital with energy Ei^^j in each dot, the Hubbard Hamiltonian can be written as 



where the operators Cjo- (q^) annihilate (create) an electron with spin a = ±1/2 in dot i. 
Further, rii^j = cj^Cia and Qi = riii + rij^ are, respectively, the spin and charge density on the 
ith dot. In Eq. ([1]), tij{B) is the matrix element describing tunneling between dots i and 
j, U is the onsite Coulomb repulsion, and V is the direct repulsion of electrons occupying 
neighboring dots. 

In the Hamiltonian ([1]) the magnetic field B is accounted for in two terms. First, it 
introduces the Zeeman splitting in the onsite energies Ei^^ = Ei + g*fiBa, with g* being the 
effective Lande factor and being the Bohr magneton. Second, it renormalizes the single- 
particle tunneling elements tij by Peierls phase factors,— such that tij{B) = tijc'^'^^'^'K For 
the three quantum dots located in the corners of an equilateral triangle we have 0i2 = 023 = 
031 = —0/3. Here, = ?)\/?)eBB? /StiTlc is the number of magnetic flux quanta threading 
the system, with e being the electron charge, c - the speed of light, % - the Planck's constant, 
and R - the distance from the center of the triangle to each dot, identified in Fig. [T](a). 

With one spin-degenerate orbital per dot we can fill the TQD with up to A^ = 6 electrons. 
To find the eigenenergies and eigenstates of A^ electrons we use the configuration interaction 
approach (CI), in which we create all possible configurations of A^ electrons on the localized 
orbitals, write the Hamiltonian Hh in a matrix form in this basis, and diagonalize the 
resulting matrix numerically.— 
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+ UY^ n^inii + -V QiQj, (1) 




i=l ij=l,i^j 
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B. LCHO-CI method 



We compare the results of the Hubbard model with a microscopic approach to the cal- 
culation of the electronic properties of a TQD starting from a confining potential, which 
we outline in this section. We start by expressing all energies in units of the effective Ryd- 
berg TZ = mle^/2e'^fi^, and lengths in the effective Bohr radius, as = eh^/mle^, where m* 
is the electron effective mass and e is the dielectric constant of the material. With GaAs 
parameters, m* = 0.067 rriQ and e = 12.4, we have TZ = 5.93 meV and = 9.79 nm. 

A single electron in the TQD in the presence of an external perpendicular magnetic field 
is described by the Hamiltonian 
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H, = (-zV + A(r))=^ + J2 V^ir) + ^^^(r) (2) 

1=1 



where A(r) is the effective vector potential, Vi{r) is the confining potential of the i-th dot, 
and VB(r) is the potential due to the additional gates, which control the potential barriers 
between dots. 

We choose the vector potential in the symmetric gauge A(r) = r2c/4(— yx + a;y) centered 
at the geometric center of the triangle of the dots.— Here the cyclotron energy = huOc/TZ 
with ujc = eB/mlc. The confining potential of each dot is approximated by a Gaussian 



(""fir) 



. Further we separate the Gaussian potential into the harmonic 



and anharmonic parts, 

V.{r) = -V,o + ^ (^)' + = Vr + SV. (3) 

where f2j^o = '^\JVifl/di is the effective characteristic energy of the harmonic confinement. In 
order to tune the height of the tunneling barrier between dots independently of the confining 
potential, we introduce Gaussian barrier potentials located between each pair of dots.— In 
the device depicted in Fig. [H^b) these potentials are generated by the gates shown as white 
regions. The coordinate system used to define the Gaussian barriers is summarized in Fig. [21 
Here we assume that the narrow gate is oriented along the axis x, which means that the 
geometry in Fig. [2] applies specifically to the lower left-hand gate of Fig. [11(b). The Gaussian 
barrier can now be defined as 

VBir) = E Vi,^\r) = vH^exp f- J- - ^) , (4) 
j=i \ ^yj / 



with the global and local coordinate systems related by 



XBj) cosT]j + {y- yBj) sinrij 



Vj = -{x - XBj) sinrjj + {y - yBj) cosrjj. (5) 

If we choose both the global gauge A(r) and a computational basis centered at the origin 
of a TQD,— i2S we find a very poor convergence of results as a function of the size of single- 
particle basis, especially for large interdot distances, when each dot should essentially be 
considered separately, with its own vector potential.— To remedy this, we divide the system 
into three regions, whose boundaries are marked in Fig. [I](a) by dashed lines, and in each 
region define the vector potential in the form 

Ai = Y [-(y - y^)± +{x- Xi)y] , (6) 

i.e., centered in the potential minimum of the respective dot. 

We solve for the eigenenergies and eigenvectors of the Hamiltonian ([2]) in the basis com- 
posed of harmonic oscillator states (HO) of each dot in the magnetic field 



(r|mm) = 0i„m(r) = exp 



-^^[-yiX + Xiy) 



(7) 



where ~ ^i) HO orbitals of ith dot, satisfying the Schrodinger equation 

:-.V + A,(r))^ + Vr] 0« (r - r,) = £,^^</.« (r - r,). (8) 
The energy associated with the HO state 0^*^, 



is defined in terms of energies VLi^± = Qi^h^^cf^, with the hybrid energy Qi^h = y ^f,o + 

The eigenfunctions of Eq. ([8]) are the Fock-Darwin (FD) orbitals, whose explicit form as a 

function of z = x + iy and z = x — iy is 

(-1)^2 ( zz\ \ (z/V2kh)"'~"' for m>n 

/27illf, V nil \2ll,J y{z/V2k,h) for n>m 



Here ni = max(n, m), n2 = min(n, m) and the hybrid length li h = yi/^i,h- Further, is 
the generalized Laguerre polynomial defined as 



The phase factor e~^^~^'^+^'^^ of the basis function (/)j„m in Eq. ([7]) is due to the gauge 
transformation Aj A. This additional phase factor depends on the magnetic field and 
the distance of dot i from the origin, and leads to the fiux-dependent factor renormalizing 
the tunneling matrix elements in the Hubbard model. 

Now we can represent the single-electron eigenvalue problem of the TQD in matrix form 
in a restricted Hilbert space formed by hq FD orbitals from each dot, with dimension Norb = 
3no, as: 

Hoa(") = e„Sa(") (12) 

where Hq is the Hamiltonian matrix for Hq in Eq. ([2]) and S is the overlap matrix due to 
the non-orthogonality of the basis and a*^") is an eigenvector corresponding to the eigenvalue 
En- The eigenstates are given by 

Norb 

Mr) = E 4"V.(r) (13) 

k=l 

where the composite index k = {inm}. The Hamiltonian and overlap matrix elements can 
be obtained efficiently if we expand the FD orbitals as linear combinations of the zero-field 
HO orbitals (pn},nyi'r) with characteristic, magnetic-field dependent energy ^1^. 

n+m 

</)„^(r) = E Ar0Km-s,s(r), (14) 

s=0 

where 



Anm _ / ^!(^ + m- S)\S\ (-1)^ 

V 2«2- ^ \ ^ , kUn - k)\(s - k)\(m - s + ky: ^ ' 

Then the integration needed to obtain the Hamiltonian and overlap matrix elements can be 
separated into x- and y-dependent parts and each integral can be carried out analytically. For 
the barrier potential Ve(r) such a separation is complicated by the appearance of an xy term 
in the exponent. This term can be eliminated by a transformation to the local coordinate 
system defined in Eq.(l5]), after which the integrals can be carried out analytically. 

The generalized eigenvalue problem formulated in Eq. f|T2|) can be cast into a standard 
eigenvalue problem 

H'b(") = e„b("\ (16) 

where H' = (yS)"iHo(yS)~^ and b(") = v^a("). The matrix is found by solving 
the eigenvalue problem SV5 = V5E5. Here V5 is the matrix of eigenvectors and E5 is 



the diagonal matrix with eigenvalues. Then ^/S is obtained by y/S = V^E^'^^Vg. The 
off-diagonal elements of the effective Hamiltonian H' correspond to the tunneling elements 
in the Hubbard model. To see how the gauge transformation automatically takes care of 
the phase change of the tunneling elements let us consider a resonant TQD system where 
all three confining potentials are identical. Since we will consider only s orbitals from each 
dot (i.e., n = m = 0), we shall use the simplified notation |iOO) = Then the off-diagonal 
matrix element of the Hamiltonian Hq are 



I^oIj) = + {^\ - ^(r - r,r\j) + {^\ E V,{r)\j) 

^ k=i 



The overlap matrix element takes the form 



(z|j) = exp 



A 4 ^ 



and the second term in Eq. (1171) is 



W - 



1j) 



02/2 



2 + 



(17) 



(18) 



(19) 



4 ^ 4 

If we neglect the three-center integrals (i|\4(r)|j) for k ^ the last term in Eq. (JTTl) is 
obtained in the form 



~2V H"^ 

(i|y,(r) + F,(r)|j) = ^^^exp 



T-xr.)-— + ^ + 
^'^^'^^ 2{2ll + d^)\:^ All^ 16 ^ 



(20) 

The common overall phase ■ (r^ x r^) is proportional to the magnetic field and the area 



of the parallelogram formed by the vectors and r^. Now Eq. 



becomes 



(21) 



where the amplitude Aij has complicated dependence on the magnetic field but is generally 
positive and decreases exponentially as the magnetic field increases. This exponential de- 
crease is due to the suppression of the overlap of orbitals from different dots, resulting from 
the decrease of the effective radius of the wave function with the increasing magnetic field. 
The off-diagonal element of the effective Hamiltonian H' differs from that in Eq. (IT7|) due 
to the existence of the overlap matrix S, but the behavior of the phase and the amplitude is 
the same. Thus the tunneling parameter in the Hubbard model in the presence of magnetic 
field acquires a field- dependent phase proportional to the flux, and amplitude which decays 
exponentially with the flux. 
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The eigenstates of the above single-electron problem are linear combinations of the har- 
monic oscillator orbitals (LCHO). We use these LCHO extended molecular orbitals to solve 
the many-electron problem of the TQD system. The Hamiltonian of this system is 

H = Y.i^^ + ef)clci^ + ^ J2 {^3\vc\kl)clc]„,ck^,ci^ (22) 

i<J ijklcra' 

where z, j, fc, / enumerate the LCHO orbitals and a, a' are spin indices. The operators 
(cjcr) create (annihilate) an electron on the spin-orbital (z,cr), while = g*m*^VLc(J /2mQ 
is the Zeeman energy. In the following discussions, the Zeeman energy is accounted for 
only in the sections corresponding to three electrons and the addition spectra, where it is 
responsible for the transition between spin polarized and spin unpolarized ground state. In 
order to make the transition more clear in the corresponding figures, we have chosen a model 
value of g* = —0.02 instead of the usual value corresponding to GaAs {g* = —0.44). 

The second term of the above Hamiltonian is scaled by Coulomb interaction matrix 
elements 

{tj\v,\ki) = jdvj dv'rd^)ri{^')y^^i,k{v'm^). (23) 

Using the Fourier transformation of the Coulomb interaction, 

{ij\v,\kl) = \j^dq fj dO^j drij*{r)e''^-'iJi{r) J drV*(r')e-^'i-^>,(r') 

1 roo r27T 

= - dq d9,{t\e'^-''\l){j\e~^'^-''\k) . (24) 

TT JO Jo 

The matrix elements of the plane wave {i\e'^'^''^\j) are evaluated analytically using the ex- 
pansions (1131 - [T^ of the LCHO orbitals in terms of the zero-field HO orbitals. The q and 
6q integrations are then carried out numerically. The Coulomb interaction matrix elements 
can be used to extract the interaction parameters in the Hubbard model. 

While the LCHO-CI approach is general, we will illustrate it on the TQD molecule with 
identical quantum dots. For a given number of electrons, we consider all possible configura- 
tions of electrons in the LCHO orbitals, calculate the Hamiltonian matrix in this configura- 
tion basis and diagonalize this matrix numerically to find the eigenstates and eigenenergies 
of the interacting many-electron TQD. In this paper, we consider single-particle basis formed 
by s orbitals from each dot and filling of the lowest electronic shell with = 1 — 6 electrons. 
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III. MAGNETIC FIELD BEHAVIOR OF THE LOWEST ELECTRONIC SHELL 



A. Magnetic field dependence of single-electron spectrum 



Let us start our analysis by discussing the single-particle spectrum of the triple dot 
molecule as a function of the magnetic field. In the basis of orbitals {|1), |2), |3)} localized on 
the respective dots, the Hubbard Hamiltonian for a single electron in the TQD on resonance 
takes a matrix form 

Ht= te'^^i^/^ E te-2-^<^/3 . (25) 
The one-electron Hamiltonian can be diagonalized by performing the Fourier transform of 

3 

the localized basis \j) into a plane wave basis \K) as \K) = J2 e^^^^~^^\j) (Ref. |3l|). The 
new basis consists of three states, with Ki = 0, K2 = 27r/3, and = —2n/3, given by: 

l^i) = 73(|l) + |2) + |3)) 



|ir2) = 45(|l) + e^'-/3|2) + e*^-/3|3) 



(26) 



_ IK3) = -1= + e-*2V3|2) + e-'^^/^'lS) 

The corresponding eigenenergies are, respectively: Ei = E — 2|t| cos (27r</)/3), E2 = E — 
2\t\ cos [27r(0 + l)/3], and E3 = E - 2\t\ cos [27r(0 - l)/3]. At zero magnetic field the three 
eigenstates form a spectrum with a non-degenerate, standing wave (zero effective angular 
momentum) K = ground state, and two degenerate excited states with K = ±27r/3 
(effective angular momentum ±1, respectively). 

In Fig. [3](a) we show these energies as a function of the flux with different lines 
corresponding to each effective angular momentum. The calculations were performed for 
model parameters E = and t = —0.01187^, and in the absence of the Zeeman energy. 
We find that the one-electron energy spectrum is composed of three levels, whose energies 
undergo Aharonov-Bohm oscillations with period = 3 flux quanta and amplitude 2\t\ 
around the single-dot energy E. At (p = {2n + l)/2 with n = 0, 1, ... we find a degenerate 
ground state and a nondegenerate excited state of the system. On the other hand, for = n 
the degeneracy is inverted, i.e., the ground state is nondegenerate while the excited state is 
doubly degenerate. The levels correspond to different quantum numbers, and hence cross 
without interaction, leading to degeneracies. 
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We tested the behavior of the energy spectrum of the Hubbard model against the micro- 
scopic LCHO approach. We assume the depth of the Gaussian potentials Vq = 5.864 71, their 
characteristic width d = 2.324 as, and the distance between dot centers |rj — rj\ = 4.85 as 
based on fitting to the electrostatic confinement produced by a model lateral gated quantum 
dot device.- As discussed in the previous section, the fiux-dependent phase factor is due to 
the gauge transformation in LCHO approach. In the inset of Fig. [3]^b) we show the single- 
particle energies as a function of the magnetic flux calculated with the LCHO method with 
only one HO orbital per dot. The resulting spectrum does exhibit the periodic degeneracies 
of levels. It differs, however, from that in Fig. [3](a) in two aspects. First, as a function of the 
magnetic field all energies undergo a diamagnetic shift towards higher energies. This shift 
is, in most part, due to the behavior of single-dot energies, which in the LCHO approach 
are ef^Q = Qh, and therefore increase with the magnetic field. In the Hubbard model, on 
the other hand, we have assumed these energies to be constant, irrespective of the number 
of flux quanta. In the main panel of Fig. [3|b) we have redrawn the LCHO spectrum with 
the diamagnetic shift removed by subtracting the reference energy 

m). (27) 

The renormalized spectrum can be directly compared with the energy spectrum of the 
Hubbard model. Both spectra oscillate with increasing magnetic field. The second difference 
between the two spectra involves the amplitude of oscillation, which remains constant in the 
Hubbard approach, but decreases in the LCHO treatment. This feature can be understood 
in terms of decrease of the magnitude of the effective tunneling parameter \t\ with increasing 
magnetic field. This can be overcome by reducing the height of tunneling barriers between 
dots using additional gates as discussed in Section HTl 

From Fig. [3] it is apparent that by adjusting the magnetic field and barrier height we 
can engineer the degeneracies of the single-particle states. This property of the triple dot 
molecule is of key importance when the system is being filled with electrons. 

B. Two electrons 

Let us start with N = 2 electrons confined in the triple dot molecule. In order to simplify 
the notation, in the following sections, unless the opposite is explicitly stated, we shall denote 
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3 

(^001 i-iV + Ai{r)f + Y.'^j{ 



the complex and magnetic field- dependent hopping parameter Uj{(j)) by tij. 

We can classify the two-electron states into singlets and triplets according to their to- 
tal spin. Let us start with the triplet subspace, with both electrons spin-down. The basis 
consists of three singly-occupied localized configurations: |Ti) = Cg^cljO), IT2) = cl^CgJO), 
IT3) = C31C2JO). Each configuration has the same energy 2E + V, and each pair of config- 
urations is coupled via the single-particle tunneling elements only. Therefore, the Hubbard 
Hamiltonian written in this basis is identical to the single electron Hamiltonian, Eq. fl25l) . 
except that all off-diagonal tunneling elements acquire a negative phase. As a result, the 
triplet eigenvectors |Ti), IT2), IT3) can be expressed as Fourier transforms of the basis states 
\Tj) in the same way the single-particle molecular orbitals are expressed in terms of lo- 
calized orbitals shown in Eq. fl2Bl) . The two electrons either move clockwise, counter- 
clockwise, or stand still. The three eigenenergies corresponding to these eigenvectors are, 
respectively, = 2E + V + 2\t\ cos (27r0/3), E^ = 2E + V + 2\t\ cos [27r(0 + l)/3], and 
E^ = 2E + V + 2\t\ cos [27r(0 — l)/3]. Note the difference in sign of t in the eigenvalues with 
respect to the single electron case. As a result, at zero magnetic field we obtain the doubly 
degenerate lowest-energy state Et = 2E + V — \t\, and a non-degenerate excited state with 
energy E^ = 2E + V + 2\t\, As the magnetic field increases, the triplet energies oscillate 
with the period of 3 flux quanta. 

Let us now move on to the singlet subspace. The singly-occupied singlet configu- 
rations 15*1), 15*2), and 15*3) are obtained from the triplet configurations |Ti), IT2), and 
IT3) by flipping the spin of one electron and properly antisymmetrizing the configura- 
tions. For example, the configuration l^i) = (c2|cj| + c||C2|) |0). In the same way, 

1*^2) = :72 ('^si'^iT ~'~ '^U'^st) I*^) 1*^3) ~ 73 ('^3J.'^2T + 4j,'^3t) III addition to the singly- 
occupied configurations there are also three doubly-occupied configurations 15*4), l^s), and 
\Sq), such that, e.g., {84) = c\^c\^\0). The singly-occupied configurations are characterized 
by energies 2E + V, and, just as the triplets, they are coupled by tunneling matrix elements. 
Unlike in the triplet case, however, these off-diagonal elements do not acquire the negative 
sign. On the other hand, the energies of all doubly-occupied configurations are 2E + U, i.e., 
contain the element describing the Coulomb onsite repulsion, making these energies larger 
than those of the singly-occupied configurations. The Hubbard Hamiltonian does not mix 
the configurations 15*4), IS^), and \Sq) with each other, but does mix the singly and doubly 
occupied subspaces. Here again it is convenient to Fourier transform the singlet basis set 
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into the form {IS"!), IS'2), jS's), |5'4), 15*5), \Sq)}. In this basis the full singlet Hamiltonian can 
be written as two 3x3 block diagonal matrix coupled through terms that account for the 
interactions between singly and doubly occupied configurations: 



JDi C 
/D2 



(28) 



where / is the 3x3 identity matrix, the vector Dj = 
[2E + V - 2\t\ cos (27r0/3) ,2E + V - 2\t\ cos (27r(0 + l)/3) ,2E + V - 2\t\ cos {2n{^ - l)/3)], 
and the vector = [2E + U,2E + U, 2E + U]. The coupling matrix C is 



cos ( ^ j gj27r(/./3 g-i27r0/3 



ig-j27r<^/3 _^ ^i2TT (2(1,-1)/ 3'^ gi27r0/3 |^]_ _^ g-47rj/3^ j27r((/.+l)/3 |^]_ _^ g~27ri/3^ 
lg-j27r(/)/3 _^ gi27r(2</,+l)/3^ gi27r</-/3 |^]^ _^ g47ri/3^ g-j27r(0-l)/3 |^]^ _^ g27ri/3^ 

In Fig. 111(a) we have plotted the low-energy spectrum for the TQD with 2 electrons, assum- 
ing the Hubbard parameters V = 0.42 TZ and U = 2.56 TZ (Ref. llTt ). The three dashed 
lines correspond to the eigenvalues of the triplet Hamiltonian while the three solid lines are 
eigenvalues of the singlet Hamiltonian. For the 6x6 singlet Hamiltonian, there is an addi- 
tional three-fold degenerate and non-oscillating eigenvalue at higher energy, originating from 
the doubly-occupied configurations (not shown in the figure). The main result, apparent 
in Fig. 111(a), is the existence of transitions between the spin singlet and triplet, occurring 
periodically as a function of the magnetic flux. In the region between = and = 1 
the ground state is a singlet for < 1/4 and > 3/4, and triplet for 1/4 < 30/4. This 
alignment of phases repeats for each subsequent flux quantum. 

Upon the inclusion of Zeeman energy we find that the intervals of stability of the singlet 
phase decrease in each subsequent period. The spin oscillations are eventually suppressed 
leading to a continuous triplet ground state at sufficiently high magnetic fields. 

The existence of spin oscillations is confirmed by results of the LCHO-CI calculation 
presented in Fig. Hl^b). The inset of Fig. Hl^b) shows that the diamagnetic shift together 
with the decrease of tunneling between dots makes it difficult to distinguish more than one 
oscillation. But after removing the diamagnetic shift by subtracting the reference energy 
Eq, defined by the ground state energy of the two-electron system without tunneling, the 
resulting energy spectrum [main panel in Fig. 111(b)] agrees well with the Hubbard model 
except for the exponential decay of the amplitude of energy oscillations. 
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C. Three electrons 



The three electron case at zero magnetic field has been analyzed in detail in Ref. Il7|. Fol- 
lowing that scheme, we start our treatment with the completely spin-polarized system, i.e., 
one with total spin S* = 3/2. In this case we can distribute the electrons on the three dots in 
only one way: one electron on each site with parallel spin, which gives a spin-polarized state 
1*23/2) = C3jC2|c|jO). This is an eigenstate of our system with energy i?3/2 = SE + 3V. Let 
us now flip the spin of one of the electrons. This electron can be placed on any orbital, and 
with each speciflc placement the remaining two spin-down electrons can be distributed in 
three ways. Altogether we can generate nine different conflgurations. Three of these config- 
urations involve single occupancy of the orbitals. They can be written as \a) = cLcLcLlO), 



\b) = c||C3|C2||0), and |c) = C2|c{|C3-|.|0). The remaining six configurations with double 



J J J 



occupancy are \d) = cl^c[^c[^\0) , \e) = cl^c[^c[^\0) , \f) = 4^4^4||0), \g) = 4^4^4110)' 
\h) 

~ 4i4i4TlO)' b) ~ 4i4i4TlO)- these configurations are characterized by the same 
projection of total spin, 5*2 = —1/2. Moreover, the doubly-occupied configurations are also 
the eigenstates of total spin, with S = 1/2, while the total spin of the singly-occupied config- 
urations is not defined. In the basis of the nine configurations we construct the Hamiltonian 
matrix by dividing the 9 configurations into three groups, each containing one of the singly- 
occupied configurations |a), \b), and |c), respectively. By labeling each group with the index 
of the spin-up electron, the Hamiltonian takes the form of a 3 x 3 matrix: 



Hl/2 



Hi 


T12 


ft 


ft 
12 


H2 


^23 


T31 


ft 

23 


^3 



(29) 



The diagonal matrix, e.g. 



Hi 



3E + 2V + U 



"-23 
''13 



^23 ~^13 

3E + 2V + U ti2 

3E + 3V 

describes the interaction of three configurations which contain spin- up electron on site 1, 
i.e., two doubly-occupied configurations \d) and |e), and a singly-occupied configuration \a). 
The remaining matrices corresponding to spin-up electrons localized on sites 2 and 3 can be 
constructed in a similar fashion. The interaction between them is given in terms of effective 
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and magnetic field dependent hopping matrix 



T- ■ 



-Uj 
-tij 
Hij 



Upon diagonalization of tlie Hamiltonian fl2^ we obtain nine levels, of which one corresponds 
to the total spin S = 3/2, and eight - to the total spin S = 1/2. The energy of the high- 
spin state is the same as that of the configuration |a3/2) discussed above, except for the 
Zeeman contribution, which is different due to the different spin projection 5*^ of the two 
configurations. 
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Let us now discuss the energy spectrum of the system at zero magnetic field. In Ref. 
we have shown that this spectrum is composed of two segments. In the low-energy region 
we find two S = 1/2, Sz = —1/2 states, which form a degenerate pair for the TQD on 
resonance. At the energy equal to 3 J/2, where J is the exchange energy, we find one 
S = 3/2, Sz = —1/2 state. These three levels are built out of singly-occupied configurations. 
The high-energy part of the spectrum consists of three pairs of states, composed of doubly- 
occupied configurations. The two parts of the spectrum are separated by an energy gap 
proportional to the onsite Coulomb element U. In the following we shall focus on the 
low-energy segment of the spectrum only, shown in Fig. [5t^a). 

Before we discuss the three-electron spectrum at finite magnetic field, we first account for 
the correct degeneracy of the energy levels by including states with all possible orientations 
of the total spin Sz- In this case we have two pairs of states with low spin: one pair with 
Sz = +1/2, and another with Sz = —1/2. These two pairs form a degenerate quadruplet at 
= 0. The high-spin state, on the other hand, is a manifold of four states, with Sz = ±3/2 
and Sz = ±1/2. Let us now consider the spectrum at finite magnetic fields accounting 
for the Zeeman energy. The quadruply degenerate low-spin state splits into two branches 
separated by the Zeeman energy, see Fig. [5]^a), reflecting the different orientations of Sz- 
Further, the energies of the states composing each pair oscillate with the magnetic field, 
and cross each other at = mr /2, n = 0, 1, . . . (not seen on this scale). These oscillations 
have a period A0 = 1, different from the period of three flux quanta present for one and 
two electrons. The amplitude is also a non-trivial function of the hopping parameters t, 
being more than two orders of magnitude smaller than As for the high-spin state, its 
four-fold degeneracy is lifted by the Zeeman energy, but the constituent levels do not exhibit 
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any oscillations. With increasing magnetic field the S = 3/2, = 3/2 spin-polarized state 
lowers its energy with respect to the ground S = 1/2, Sz = 1/2 state, and at a critical value 
of the magnetic field becomes the ground state. The critical magnetic field is given by the 
condition g*^Bc = 3 J/2, under which the Zeeman energy equals the exchange energy. The 
results of the Hubbard calculations are in agreement with the three-electron energy spectra 
calculated within the LCHO-CI approach, shown in Fig. [5]^b). Note that a similar analysis 
was reported in Ref. for the TQD composed of shallower, parabolic dots with larger 
interdot tunneling. These calculations revealed an additional total spin oscillation between 
the = 1/2 and S = 3/2 phases, which occurs for the same Sz component. We were able to 
reproduce this oscillation within the LCHO-CI approach using shallower dots, but not within 
the Hubbard model. This is because the spin oscillation is due to the magnetically-induced 
reduction of the interdot tunneling element. However, the importance of these oscillations 
is minor due to the dominant role of the Zeeman energy. As discussed above, the Zeeman 
energy leads to the onset of a spin polarized phase at a sufficiently high magnetic field, which 
suppresses any spin oscillations. 

D. Four electrons 

The four-electron configurations correspond to two holes, created when two electrons are 
removed from the fiUed-shell configuration. With two holes we can form only the spin singlet 
and triplet configurations of the system. 

Let us focus on the triplets first. They involve one electron spin- up occupying the first, 
second, or third dot in the presence of an inert core of three spin-down electrons. If we 
denote the operator creating a hole on dot i with spin a by h]^, we can write the three 
basis configurations in this subspace in the form {T^;^^) = h\^hl^\Ne = 6) = CgiCgj^Caj^c^jO), 
\T2^^) = hl^h\^\Ne = 6) and \T^^^) = hl^hl^\Ne = 6). The two-hole triplet Hamiltonian is 
given by 

Ue + U + 5V t23 hs 
Ht= t*3 4E + U + 5V ti2 ■ (30) 

i*3 t*2 AE + U + W 

Note that the above Hamiltonian differs from that describing the two-electron triplet 
subspace in that the off-diagonal tunneling matrix elements do not acquire the additional 
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negative phase. 

Let us move on to the two-hole singlet configurations. The singly-occupied states involve 
the two holes occupying two different dots, while the doubly-occupied states hold both holes 
on the same dot. This situation is analogous to the two-electron case described earlier, and 
the two-hole singlet Hamiltonian written in the appropriately rotated basis is analogous to 
that shown in Eq. (!28l) . The only difference is that the diagonal vectors Di and D2 will 
contain two-hole, instead of two-electron energies: AE + U + 5V instead of 2E + V for 
singly-occupied configurations, and + 2U + W instead of 2E + U for doubly-occupied 
configurations. Also, the factors |t| in Di acquire the opposite sign, while this additional 
phase does not appear in the coupling matrix C for the two-hole case. 



In Ref. [171 we have diagonalized the singlet and triplet Hamiltonians at zero magnetic 
field. We found that in this case the total spin of the two-hole ground state depended on 
the interplay of Hubbard parameters. For a typical case of 2|t| < U — V the ground state is 
a spin triplet. Clearly, the appearance of the finite magnetic moment of the ground state is 
made possible by the degeneracy of the single-particle excited state at zero magnetic field. 
As we increase the field this degeneracy is removed, so we may expect a transition to a spin 
singlet. The periodic reappearance of the degeneracy should lead to spin oscillations. This 
is indeed what we observe in the energy spectrum, whose low-energy segment is plotted in 
Fig. int^a) as a function of the number of magnetic fiux quanta. Again, singlet eigenvalues 
have been plotted with solid lines and dashed lines for the triplets. As for the case of = 2, 
transitions between triplet and singlet ground state appear as we increase the magnetic fiux, 
except that the alignment of phases seen for N = 2 electrons is inverted. Furthermore, as in 
the case of two electrons, the introduction of the Zeeman term will favor the spin alignment 
of the triplet configuration, suppressing the triplet-singlet transitions at high magnetic field. 
These predictions of the Hubbard model are confirmed by the LCHO-CI calculation, whose 
results are presented in Fig. [6](b). Again, the original spectrum is shown in the inset, while 
the main panel shows the energies without the diamagnetic shift. 

E. Five electrons 

Five electrons correspond to a single hole. The single-hole Hamiltonian can be obtained 
from the single-electron Hamiltonian by appropriately modifying the diagonal terms and 



17 



setting tij ^ —tij (Ref. llTI ). For the triangular triple dot on resonance this symmetry is 
reflected in the energy spectrum of the hole as shown in FiglSt^a). For the one- hole problem at 
zero magnetic field, the opposite sign of the off-diagonal element leads to a doubly- degenerate 
hole ground state. This behavior is confirmed in the LCHO-CI calculations, whose results 
are shown in Fig. MJo)- 



IV. CHARGING DIAGRAM OF THE RESONANT TRIPLE DOT 

We can now construct the charging diagram of the triple dot molecule as a function of the 
magnetic field. For any number of electrons (1 to 6) and any quantum dot energy E, we 
obtain the ground-state energy Egs{N) and the corresponding total spin by diagonalizing 
the Hubbard Hamiltonian. We use these energies to calculate the chemical potential of the 
triple quantum dot molecule fi{N) = EcsiN + 1) — Egs{N). When fi{N) equals the chem- 
ical potential /xl of the leads, the N + 1st electron is added to the A^-electron quantum-dot 
molecule. This establishes the total number of electrons in the quantum dot molecule and 
their total spin as a function of quantum-dot energy E relative to the chemical potential of 
the leads fii- In the case of LCHO-CI approach, the relevant quantum dot energy is the 
single-particle reference energy E. Changes in electron numbers can be detected by Coulomb 
blockade (CB), spin blockade, or charging spectroscopies.-^i^ The calculated stability dia- 
gram, with Hubbard parameters used in the previous section and taking into account the 
contribution of the Zeeman term, is shown in Fig. [8](a), while Fig. [8](b) shows the stability 
diagram computed using the LCHO-CI approach. Note that in this figure the oscillations 
of the stability lines corresponding to the condition fi{N) — /i^ = are not visible due to 
the energy scale. As explained in the previous sections, the differences among both addition 
spectra are due to the diamagnetic shift and the suppression of the inter-dot tunneling. 
These effects appear naturally in the LCHO-CI approach but are not taken into account in 
the Hubbard model. 

Let us explain the addition spectrum as we change the single-dot energy E of each dot 
with respect to the chemical potential of the leads fii = 0. From the condition fii = E — 2\t\, 
the energy Egs{^) corresponding to the addition of the first electron can be approximated 
by E{1) ^ 2\t\. At this energy the first Coulomb blockade peak of the triple quantum 
dot molecule should be observed. Similar arguments based on the results of the previous 
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section can be applied to find the CB peaks corresponding to the addition of the remaining 
electrons. Note that the prominent energy gap that appears for the addition of the fourth 
electron on the TQD is due to the large on-site Coulomb repulsion U. The ground state 
of the three-electron TQD corresponds to one electron occupying each dot, and therefore 
the addition of a new electron will increase the energy of the system by order of U. This 
term does not appear in any other addition processes, in which only the interdot Coulomb 
element V is relevant. 

The spin oscillations in the system with two and four electrons, as well as the spin 
transition in the three electron system, lead to a strong modulation of the current through 
the TQD. An schematic representation of the amplitude of Coulomb blockade peaks as a 
function of the magnetic flux is given in Fig. [HI for different numbers of electrons confined 
in the triple dot. Here we assume that the leads are spin-unpolarized, and the transport 
involves only the lowest-energy level of the molecule. The vertical arrow shows the transition 
in the three electron system from S = 1 /2 to S = 3/2 state. When the two electrons are in 
a spin triplet with S = 1, adding an electron can create a spin polarized final S = 3/2 three 
electron droplet, and the current is high. However, if the two-electron system is in a singlet 
5 = state, the final S = 3/2 state cannot be reached by adding a single electron, and 
the current is spin blockaded. Hence quenching of the tunneling current in the two-electron 
droplet is a signature of a spin-polarized three-electron state and a spin-singlet two-electron 
state. Adding a fourth electron to a three-electron droplet is equivalent to adding a hole to 
a two-hole droplet. Hence the oscillation in CB peak amplitude, but shifted in phase since 
the holes start as triplets and electrons start as singlets. 



V. CONCLUSIONS 



In conclusion, we have studied the effect of the magnetic field on the electronic properties 
of a triple triangular quantum dot molecule. Exact results for few-electron spectra in a 
magnetic field where obtained for identical dots in the Hubbard model. Aharonov-Bohm 
oscillations of a single electron, singlet-triplet spin oscillations for pairs of electrons and 
pairs of holes, as well as a transition from the frustrated magnetic state to spin polarized 
state for a half-filled lowest electronic shell are predicted. The impact of spin transitions 
on the stability diagram and modulation of the current through the TQD molecule with 
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increasing magnetic flux are discussed. The results of the Hubbard model are supported by 
a microscopic calculation using a general LCHO-CI approach extended to finite magnetic 
fields. 
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FIG. 1: (a) Schematic image of the triple-dot molecule. Symbols denote the Hubbard parameters 
for the system on resonance. The magnetic field is perpendicular to the plane of the molecule, 
(b) Contour plot of the potential created in the lateral triple-dot device by a typical layout of the 
metallic gates. White zones show schematically the gates controlling the barriers between dots. 
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FIG. 2: Schematic picture of the potential created by one of the Gaussian barrier centered at 
{xBj,yBj) and the relation between the original variables {x,y) and the system {xj,yj). 
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FIG. 3: Single-particle energies of the triple dot molecule as a function of the number of magnetic 
flux quanta calculated within the Hubbard model (a), and with the LCHO approximation (b). In 
figure (b) the inset corresponds to the actual spectrum while in the main panel the diamagnetic 
shift has been removed to facilitate the comparison with results from Hubbard model. 
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FIG. 4: (a) Energy spectrum of the TQD with A'' = 2 as a function of the number of magnetic 
flux quanta calculated within the Hubbard model. Solid lines correspond to the triplet levels, 
while dashed lines correspond to the singlets, (b) The same energy spectrum calculated with the 
LCHO-CI approximation after removing the diamagnetic shift. Inset corresponds to the actual 
spectrum with diamagnetic shift. 
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FIG. 5: (a) Energy spectrum of the TQD with three electrons versus the dimensionless flux ^ 
including the Zeeman term. All energies are measured from the ground state energy at = 0. (b) 
The same spectrum calculated using the LCHO-CI technique without diamagnetic shift. 
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FIG. 6: (a) Energy spectrum of the TQD with four electrons versus the dimensionless flux (p in the 
absence of Zeeman energy. Triplet levels plotted with dashed lines while singlets are indicated with 
solid lines, (b) The four-electron spectrum obtained using the LCHO-CI technique with (inset) 
and without (main panel) diamagnetic shift. 
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FIG. 7: (a) Energy spectrum of the TQD with five electrons versus the dimensionless flux 0, 
plotted with the zero-field ground-state energy treated as reference, (b) The five-electron spectrum 
obtained using the LCHO-CI method with (inset) and without (main panel) diamagnetic shift. 
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FIG. 8: (Color online) (a) Charging diagram for the TQD versus the number of magnetic flux 
quanta (p calculated in the Hubbard model. Colors indicate the total spin of the ground state: 
S = (blue), 5 = 1/2 (cyan), 5=1 (yellow), and 5 = 3/2 (red), (b) The diagram obtained using 
the LCHO-CI method. 
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FIG. 9: Qualitative scheme of the current passing through the TQD filled with A'" electrons ex- 
tracted from the Hubbard results. Thick solid lines correspond to a large tunneling current while 
thin-dotted lines correspond to a small current under the spin blockade condition. 
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